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Recent theories and experiments have suggested that strong spin— orbit cou- 
pling effects in certain band insulators can give rise to a new phase of quan- 
tum matter, the so-called topological insulator, which can show macroscopic 
entanglement effects [1, 2, 3, Q, 51,1a, 2]. Such systems feature two-dimensional 
surface states whose electrodynamic properties are described not by the con- 
ventional Maxwell equations but rather by an attached axion field, origi nally 
proposed to describe strongly interacting particles 

It has been proposed that a topological insulator [2j with a single spin-textured 
Dirac cone interfaced with a superconductor can form the most elementary unit 
for performing fault-tolerant quantum computation jl^ . Here we present an 
angle-resolved photoemission spectroscopy study and first-principle theoretical 
calculation-predictions that reveal the first observation of such a topological 
state of matter featuring a single-surface-Dirac-cone realized in the naturally 
occurring Bi2Se3 class of materials. Our results, supported by our theoreti- 
cal predictions and calculations, demonstrate that undoped compound of this 
class of materials can serve as the parent matrix compound for the long-sought 
topological device where in-plane surface carrier transport would have a purely 
quantum topological origin. Our study further suggests that the undoped com- 
pound reached via n-to-p doping should show topological transport phenomena 
even at room temperature. 



Note added: The counter- doping method to reach the fully undoped compound - the 
bulk-insulator state, proposed in this paper, has been achieved in a recent work in Nature: 
[http: / /dx.doi.org/10. 1038 /nature08234| (2009)] 
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It has been experimentally shown that spin-orbit coupling can lead to new phases of 
quantum matter with highly nontrivial collective quantum effects (4, 5, 6]. Two such phases 
are the quantum spin Hall insulator 4] and the strong topological insulator 0, [(], Ql both 
realized in the vicinity of a Dirac point but yet quite distinct from graphene [16] . The 
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ator phase contains surface states (SSs) with novel electromagnetic 
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151 ] . It is currently believed that the Bii^Sbr insu- 



lating alloys realize the only known topological-insulator phase in the vicinity of a three- 



dimensional Dirac point |5J , which can in principle be used to study topological electromag- 



netic and interface superconducting properties [§J, |9j, llO, ll4j . However, a particular challenge 
for the topological-insulator Bii^Shj; system is that the bulk gap is small and the material 
contains alloying disorder, which makes it difficult to gate for the manipulation and control 
of charge carriers to realize a device. The topological insulator Bii^Sbr features five surface 
bands, of which only one carries the topological quantum number 6J. Therefore, there is 
an extensive world-wide search for topological phases in stoichiometric materials with no 
alloying disorder, with a larger gap and with fewer yet still odd-numbered SSs that may 
work as a matrix material to observe a variety of topological quantum phenomena. 

The topological-insulator character of Bix^Sbz [5j, |6| led us to investigate the alternative 
Bi-based compounds Bi 2 X 3 (X=Se, Te). The undoped Bi 2 Se3 is a semiconductor that be- 
longs to the class of thermoelectric materials Bi 2 X 3 with a rhombohedral crystal structure 
(space group D 5 3 d (R3m); refs 17,18). The unit cell contains five atoms, with quintuple lay- 
ers ordered in the Se(l)-Bi-Se(2)-Bi-Se(l) sequence. Electrical measurements report that, 
although the bulk of the material is a moderately large-gap semiconductor, its charge trans- 
port properties can vary significantly dep ending on the sample preparation conditions 19) , 



with a strong tendency to be n-type 



20, 



2l| owing to atomic vacancies or excess selenium. 



An intrinsic bandgap of approximately 0.35 eV is typically measured in experiments [22j, |23|, 
whereas theoretical calculations estimate the gap to be in the range of 0.24-0.3 eV (refs 20, 
24). 

It has been shown that spin-orbit coupling can lead to topological effects in materials 

nnriQ 

that determine their spin Hall transport behaviors |4J, [5J, [6|, I7j] . Topological quantum proper- 
ties are directly probed from the nature of the electronic states on the surface by studying 
the way surface bands connect the material's bulk valence and conduction bands in mo- 
mentum space [5, [sl, 3]. The surface electron behavior is intimately tied to the number 
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FIG. 1: Theoretical calculations and experimental results : Strong spin-orbit interac- 
tion gives rise to a single-surface-state-Dirac-cone on the surface of the topological 
insulator. 




FIG. 2: Transverse- momentum k z dependence of topological Dirac bands near V. 

of bulk band inversions that exist in the band structure of a material [7]. The origin of 
topological Z2 order in Bii„ a; Sb a ; is bulk-band inversions at three equivalent L-points 
whereas in Bi2Se3 only one band is expected to be inverted, making it similar to the case in 
the two-dimensional quantum spin Hall insulator phase. Therefore, a much simpler surface 
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FIG. 3: The topology of the surface Dirac cone Fermi surface and the associated Berry's 
phase. 

spectrum is naturally expected in Bi 2 Se 3 . All previous experimental studies of Bi 2 Se 3 have 
focused on the material's bulk properties; nothing is known about its SSs. It is this key ex- 
perimental information that we provide here that, for the first time, enables us to determine 
its topological quantum class. 

The bulk crystal symmetry of Bi 2 Se3 fixes a hexagonal Brillouin zone (BZ) for its (111) 
surface (Fig. Id) on which M and T are the time-reversal invariant momenta (TRIMs) 
or the surface Kramers points. We carried out high-momentum-resolution angle-resolved 
photoemission spectroscopy (ARPES) measurements on the (111) plane of naturally grown 
Bi 2 Se 3 (see the Methods section). The electronic spectral weight distributions observed 
near the T point are presented in Fig. la-c. Within a narrow binding-energy window, a 
clear V-shaped band pair is observed to approach the Fermi level (E_p). Its dispersion or 
intensity had no measurable time dependence within the duration of the experiment. The 
'V bands cross Ej? at 0.09 A~ l along T-M and at 0.10 A -1 along T-K, and have nearly equal 
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band velocities, approximately 5 x 10 5 ms _1 , along the two directions. A continuum- like 
manifold of states - a filled U-shaped feature - is observed inside the V-shaped band pair. 
All of these experimentally observed features can be identified, to first order, by a direct 
one-to-one comparison with the LDA band calculations. Figure If shows the theoretically 
calculated (see the Methods section) (lll)-surface electronic structure of bulk Bi 2 Se 3 along 
the K — f — M /c-space cut. The calculated band structure with and without SOC are 
overlaid together for comparison. The bulk band projection continuum on the (111) surface 
is represented by the shaded areas, blue with SOC and green without SOC. In the bulk, 
time-reversal symmetry demands E(k, |) = E(—k, j) whereas space inversion symmetry 
demands E(k, |) = E(—k, j). Therefore, all the bulk bands are doubly degenerate. However, 
because space inversion symmetry is broken at the terminated surface in the experiment, SSs 
are generally spin-split on the surface by spin-orbit interactions except at particular high- 
symmetry points-the Kramers points on the surface BZ. In our calculations, the SSs (red 



dotted lines) are doubly degenerate only at T (Fig. 1 



spin-orbit-coupled material surfaces such as gold [25|, |26|] or Bii-^Sb^ (ref. 5). In Bi 2 Se3, 



'his is generally true for all known 



the SSs emerge from the bulk continuum, cross each other at f , pass through the Fermi 
level {Ep) and eventually merge with the bulk conduction-band continuum, ensuring that 
at least one continuous band-thread traverses the bulk bandgap between a pair of Kramers 
points. Our calculated result shows that no surface band crosses the Fermi level if SOC is 
not included in the calculation, and only with the inclusion of the realistic values of SOC 
(based on atomic Bi) does the calculated spectrum show singly degenerate gapless surface 
bands that are guaranteed to cross the Fermi level. The calculated band topology with 
realistic SOC leads to a single ring-like surface FS, which is singly degenerate so long as the 
chemical potential is inside the bulk bandgap. This topology is consistent with the Z 2 = — 1 
class in the Fu-Kane-Mele classification scheme ~\. 

A global agreement between the experimental band structure (Fig. la-c) and our the- 
oretical calculation (Fig. If) is obtained by considering a rigid shift of the chemical po- 
tential by about 200 meV with respect to our calculated band structure (Fig. If) of the 
formula compound Bi 2 Se3. The experimental sign of this rigid shift (the raised chemical 
potential) corresponds to an electron doping of the Bi 2 Se3 insulating formula matrix (see 
Supplementary Information). This is consistent with the fact that naturally grown Bi 2 Se3 
semiconductor used in our experiment is n-type, as independently confirmed by our trans- 
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port measurements. The natural doping of this material, in fact, comes as an advantage in 
determining the topological class of the corresponding undoped insulator matrix, because 
we would like to image the SSs not only below the Fermi level but also above it, to exam- 
ine the way surface bands connect to the bulk conduction band across the gap. A unique 
determination of the surface band topology of purely insulating Bix^Sb^ (refs 5, 6) was 
clarified only on doping with a foreign element, Te. In our experimental data on Bi2Se3, 
we observe a V-shaped pure SS band to be dispersing towards Ep, which is in good agree- 
ment with our calculations. More remarkably, the experimental band velocities are also 
close to our calculated values. By comparison with calculations combined with a general 
set of arguments presented above, this V-shaped band is singly degenerate. Inside this 'V 
band, an electron-pocket-like U-shaped continuum is observed to be present near the Fermi 
level. This filled U-shaped broad feature is in close correspondence to the bottom part of 
the calculated conduction band continuum (Fig. If). Considering the n-type character of 
the naturally occurring Bi 2 Se3 and by correspondence to our band calculation, we assign the 
broad feature to correspond roughly to the bottom of the conduction band. 

To systematically investigate the nature of all the band features imaged in our data, 
we have carried out a detailed photon-energy dependence study, of which selected data 
sets are presented in Fig. 2a,b. A modulation of incident photon energy enables us to 
probe the k 2 dependence of the bands sampled in an ARPES study (Fig. 2c), allowing for 
a way to distinguish surface from bulk contributions to a particular photoemission signal 
0. Our photon-energy study did not indicate a strong k z dispersion of the lowestlying 
energy bands on the 'U', although the full continuum does have some dispersion (Fig. 2). 
Some variation of the quasiparticle intensity near E^ is, however, observed owing to the 
variation of the electron-photon matrix element. In light of the k^-dependence study (Fig. 
2b), if the features above -0.15 eV were purely due to the bulk, we would expect to observe 
dispersion as k z moved away from the T-point. The lack of strong dispersion yet close 
one-to-one correspondence to the calculated bulk band structure suggests that the inner 
electron pocket continuum features are probably a mixture of surface-projected conduction- 
band states, which also includes some band-bending effects near the surface and the full 
continuum of bulk conduction-band states sampled from a few layers beneath the surface. 
Similar behavior is also observed in the ARPES study of other semiconductors [27]. In 
our kz-dependent study of the bands (Fig. 2b) we also observe two bands dispersing in k z 
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that have energies below -0.3 eV (blue dotted bands), reflecting the bulk valence bands, 
in addition to two other non- dispersive features associated with the two sides of the pure 
SS Dirac bands. The red curve is measured right at the T-point, which suggests that the 
Dirac point lies inside the bulk bandgap. Taking the bottom of the 'IT band as the bulk 
conduction-band minimum, we estimate that a bandgap of about 0.3 eV is realized in the 
bulk of the undoped material. Our ARPES estimated_bandgap is in good agreement with 
the value deduced from bulk physical measurements 



2311 and from other calculations that 



report the bulk band structure [20|, |24[ . This suggests that the magnitude of band bending 



near the surface is not larger than 0.05 eV. We note that in purely insulating Bi2Se3 the 
Fermi level should lie deep inside the bandgap and only pure surface bands will contribute 
to surface conduction. Therefore, in determining the topological character of the insulating 
Bi2Se3 matrix the 'IT feature is not relevant. 

We therefore focus on the pure SS part. The complete surface FS map is presented in 
Fig. 3. Figure 3a presents electron distribution data over the entire two-dimensional (111) 
surface BZ. All the observed features are centered around T. None of the three TRIMs 
located at M are enclosed by any FS, in contrast to what is observed in Bii^Sb^. (ref. 
5). The detailed spectral behavior around T is shown in Fig. 3b, which was obtained with 
high momentum resolution. A ring-like feature formed by the outer 'V pure SS band (a 
horizontal cross-section of the upper Dirac cone in Fig. 1) surrounds the conduction-band 
continuum centered at f. This ring is singly degenerate from its one-to-one correspondence 
to band calculation. An electron encircling the surface FS that encloses a TRIM or a Kramers 
ooint obtains a geometrical quantum phase (Berry phase) of n mod 2tt in its wavefunction 

7|. Therefore, if the chemical potential (Fermi level) lies inside the bandgap, as it should 
in purely insulating Bi2Se3, its surface must carry a global n mod 27r Berry phase. In 
most spin-orbit materials, such as gold (Auflll]), it is known that the surface FS consists 
of two spin-orbit-split rings generated by two singly degenerate parabolic (not Dirac-like) 

Dands that are shifted in momentum space from each other, with both enclosing the f-point 
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26] . The resulting FS topology leads to a 2tt or Berry phase because the phases from 



the two rings add or cancel. This makes gold-like SSs topologically trivial despite their 
spin-orbit origin. 

Our theoretical calculation supported by our experimental data suggests that in insulating 
Bi2Se3 there exists a singly degenerate surface FS which encloses only one Kramers point 
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on the surface Brillouin zone. This provides evidence that insulating Bi2Se3 belongs to 
the Z2 = — 1 topological class in the Fu-Kane-Mele topological classification scheme for 
band insulators. On the basis of our ARPES data we suggest that it should be possible 
to obtain the fully undoped compound by chemically hole-doping the naturally occurring 
Bi 2 Se3, thereby shifting the chemical potential to lie inside the bulk bandgap. The surface 
transport of Bi 2 Se3 prepared as such would therefore be dominated by topological effects as 
it possesses only one Dirac fermion that carries the non-trivial Z 2 index. The existence of a 
large bulk bandgap (0.3 eV) within which the observed Z 2 Dirac fermion state lies suggests 
the realistic possibility for the observation of topological effects even at room temperature 
in this material class. Because of the simplest possible topological surface spectrum realized 
in Bi 2 Se3, it can be considered as the 'hydrogen atom' of strong topological insulators. Its 
simplest topological surface spectrum would make it possible to observe and study many 
exotic quantum phenomena predicted in topological field theories, such as the Major ana 
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fermions [I4j . magnetic monopole image [9j, llOj or topological exciton condensates [15|], by 
transport probes. 

NOTE ADDED: The counter-doping method to reach the fully undoped 
compound - the bulk-insulator state, proposed in this paper, and the spin- 
texture measurements have been achieved in a recent work in Nature: 

[http: / /dx.doi.org/10. 1038 /nature08234| (2009)] 



Methods 



Theoretical calculations. The theoretical band calculations were performed with the 
LAPW method in slab geometry using the WIEN2K package [281 ] . The generalized gradient 
approximation of Perdew, Burke and Ernzerhof [29J] was used to describe the exchange- 
correlation potential. SOC was included as a second variational step using scalar-relativistic 
eigenfunctions as basis after the initial calculation was converged to self-consistency. The 
surface was simulated by placing a slab of 12 quintuple layers in vacuum. A grid of 
21x21x1 points was used in the calculations, equivalent to 48 k points in the irreducible 
BZ and 300 k points in the first BZ. To calculate the k z of the ARPES measurements 
(k z = (1 / hbar) ^j2m(Ei : i ri cos 2 9 + Vo)), an inner potential Vo of approximately 11.7 eV was 
used, given by a fit on the ARPES data at normal emission. 
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Experimental methods. Single crystals of Bi2Se3 were grown by melting stoichio- 
metric mixtures of high-purity elemental Bi and Se in a 4-mm-inner-diameter quartz tube. 
The sample was cooled over a period of two days, from 850 to 650 °C, and then annealed 
at this temperature for a week. Single crystals were obtained and could be easily cleaved 
from the boule. High-resolution ARPES measurements were then performed using 17-45 
eV photons on beamline 12.0.1 of the Advanced Light Source at the Lawrence Berkeley Na- 
tional Laboratory and beamline 5-4 at the Stanford Synchrotron Radiation Laboratory. The 
energy and momentum resolutions were 15 meV and 1.5% of the surface BZ respectively 
using a Scienta analyser. The samples were cleaved in situ between 10 and 55 K under 
pressures of less than 5xl0~ n torr, resulting in shiny flat surfaces. The surface band quasi- 
particle signal is stable throughout the entire measurement duration. The counter-doping 
method to reach the fully undoped compound - the bulk-insulator state, proposed in this 
paper, and the spin-texture measurements have been achieved in a recent work in Nature: 
|http: / /dx.doi.org/10. 1038 /nature08234| (2009)] 
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FIG. 4: [Enlarged version of Fig-1] Theoretical calculations and experimental results : Strong 
spin-orbit interaction gives rise to a single-surface-state-Dirac-cone on the surface of the topological 
insulator. 
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FIG. 6: [Enlarged version of Fig-3] The Z2 topology of the surface Dirac cone Fermi surface and 
the associated Berry's phase. 
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Fig. 1 (Caption). Strong spin-orbit interaction gives rise to a single SS Dirac 
cone. Theory (see the Methods section for calculation methods) versus experiments. a,b, 
High- resolution ARPES measurements of surface electronic band dispersion on Bi 2 Se 3 (111). 
Electron dispersion data measured with an incident photon energy of 22 eV near the F-point 
along the F — M (a) and F — K (b) momentum-space cuts, c, The momentum distribution 
curves corresponding to a suggest that two surface bands converge into a single Dirac point 
at F. The V-shaped pure SS band pair observed in a-c is nearly isotropic in the momentum 
plane, forming a Dirac cone in the energy- /c-^ ky space (where k x and k y are in the Y — K 
and F — M directions, respectively). The U-shaped broad continuum feature inside the V- 
shaped SS corresponds roughly to the bottom of the conduction band (see the text), d, A 
schematic diagram of the full bulk three-dimensional BZ of Bi 2 Se3 and the two-dimensional 
BZ of the projected (111) surface, e, The surface Fermi surface (FS) of the two-dimensional 
SSs along the K — F — M momentum-space cut is a single ring centered at F if the chemical 
potential is inside the bulk bandgap. The band responsible for this ring is singly degenerate 
in theory. The TRIMs on the (111) surface BZ are located at F and the three M points. The 
TRIMs are marked by the red dots. In the presence of strong spinorbit coupling (SOC), 
the surface band crosses the Fermi level only once between two TRIMs, namely F and 
M; this ensures the existence of a ir Berry phase on the surface, f, The corresponding 
local density approximation (LDA) band structure (see the Methods section). Bulk band 
projections are represented by the shaded areas. The band-structure topology calculated in 
the presence of SOC is presented in blue and that without SOC is in green. No pure surface 
band is observed to lie within the insulating gap in the absence of SOC (black lines) in the 
theoretical calculation. One pure gapless surface band is observed between F and M when 
SOC is included (red dotted lines). 
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Fig. 2 (Caption). Transverse-momentum k z dependence of topological Dirac 
bands near f . a, The energy dispersion data along the f-M cut, measured with the 
photon energy of 21 eV (corresponding to 0.3 /c-space length along T-Z || k z ), 19 eV (r) 
and 31 eV (-0.4 /c-space length along T-Z of the bulk three-dimensional T BZ) are shown. 
Although the bands below -0.4 eV binding energy show strong k z dependence, the linearly 
dispersive Dirac-like bands and the U-shaped broad feature show weaker k z dispersion. The 
Dirac point is observed to lie inside the bulk bandgap. A careful look at the individual 
curves reveals some k z dependence of the U-shaped continuum (see b for details), b, The 
energy distribution curves obtained from the normal-emission spectra measured using 15-31 
eV photon energies reveal two dispersive bulk bands below -0.3 eV (blue dotted lines). This 
is in addition to the two non- dispersive peaks from the Dirac-cone bands inside the gap. 
The Dirac band intensity is strongly modulated by the photon energy changes due to the 
matrix-element effects (which is also observed in BiSb; ref. 5). c, A k-space map of locations 
in the bulk three-dimensional BZ scanned by the detector at different photon energies over 
a theta (9) range of ±30°. This map (k z , k y , E v h oton ) was used to explore the k z dependence 
of the observed bands. 
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Fig. 3 (Caption). The topology of the surface Dirac cone Fermi surface and 
the associated Berry's phase, a, The observed surface FS of Bi 2 Se 3 consists of a small 
electron pocket around the center of the BZ, f. b, High-momentum-resolution data around 
r reveal a single ring formed by the pure SS V-shaped Dirac band. For the naturally 
occurring Bi 2 Se3, the spectral intensity in the middle of the ring is due to the presence of 
the 'IT feature, which roughly images the bottom of the conduction-band continuum (see 
the text). The observed topology of the pure surface FS of Bi 2 Se 3 is different from that of 
most other spin-orbit materials such as gold (Au(lll)). c, The Au(lll) surface FS features 
two rings (each non-degenerate) surrounding the f point. An electron encircling the gold 
FS carries a Berry phase of zero, characteristic of a trivial band insulator or metal, and can 
be classified by Z 2 = +1 (ref. 7). The single surface FS observed in Bi 2 Se3 is topologically 
distinct from that of gold. The single non-degenerate surface FS enclosing a Kramers point 
(f) constitutes the key signature of a topological-insulator phase characterized by Z 2 = — 1. 
d, e, Schematic SS topologies in gold and Bi 2 Se 3 for direct comparison. In gold, the chemical 
potential can be continuously tuned to be placed inside the interband gap, making the SSs 
fully gapped. In Bi 2 Se3, however, the Dirac structure of the SSs required by the Kramers 
degeneracy and time-reversal invariance ensures that they remain gapless independent of the 
location of the chemical potential within the bulk gap. If the chemical potential is placed 
inside the gap, as it would naturally lie in the purely insulating undoped Bi 2 Se 3 , the surface 
transport would be dominated by the single Dirac fermion, which is of purely topological 
origin. 
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